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JET POWER EXTRACTED FROM ADAF AND THE APPLICATIONS TO X-RAY BINARIES AND RADIO 

GALAXY FRI/FR II DICHOTOMY 

QlNGWEN Wu' AND XiNWU CAO^'^ 
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ABSTRACT 

We calculate the jet power of the classical Blandford-Znajek (BZ) model and hybrid model developed by Meier 
based on the global solutions of advection dominated accretion flows (ADAFs) surrounding Kerr black holes. 
We find that the jet power of the hybrid model is larger than that of the pure BZ model. The jet power will 
dominate over the accretion power, and the objects will enter into "jet-power-dominated advective systems," 
when the accretion rate is less than a critical value m^ = M^/MEAd, where 3 x 10"^ < m^ ^ 5 x lO""* is a function 
of black hole spin parameter. The accretion power will be dominant when m > m^ and the objects will enter 
into "flccref/on-power-dominated advective systems." This is roughly consistent with that constrained from the 
low/hard-state black hole X-ray binaries (e.g.. Fender et al.). 

We calculate the maximal jet power as a function of black hole mass with the hybrid jet formation model, and 
find it can roughly reproduce the dividing line of the Ledlow-Owen relation for FR I/FR II dichotomy in the jet 
power-black hole(BH) mass plane (2jet-A/BH) if the dimensionless accretion rate m ^ 0.01 and BH spin parameter 
j ^ 0.9-0.99 are adopted. This accretion rate m ^ 0.01 is consistent with that of the critical accretion rate for the 
accretion mode transition of a standard disk to an ADAF constrained from the state transition of X-ray binaries. 
Our results imply that most FR I galaxies may be in the ADAF accretion mode similar to the low/hard- state XRBs. 

Subject headings: accretion, accretion disks-black hole physics-galaxies: jets-X-rays: binaries-MHD 



1. INTRODUCTION 

Black hole accretion is thought to power active galactic nu- 
clei (AGNs) and X-ray binaries (XRBs). Both the UV/optical 
bumps observed in luminous quasars and the soft X-rays ob- 
served in the high/soft state XRBs can be naturally interpreted 
as blackbody emission from a cold, optically thick, and geomet- 
rically thin standard disk (SSD; e.g., Shakura & Sunvaev 1973t 
I Sun & Malkanll989h . A hot, optically thin, geometrically thick 
advection dominated accretion flow model has been developed 
in the last several decades (ADAF, or "radiative inefficient ac- 
cretion flows"; e.g., Ichimaru 1977; Rees et al. 1982; Narayan 
& Yi 1994, 1995; Abramowicz et al. 1995; see Kato et al. 2008 
and Narayan & McClintock 2008 for reviews), which can suc- 
cessfully explain most features of the nearby low-luminosity 
AGNs and low/hard-state XRBs (see Remillard & McCHntock 
2006; Done et al. 2007; Yuan 2007; Ho 2008 for recent re- 
views). Their spectral energy distributions (SEDs) can be well 
reproduced by the ADAF model. It was found that the hard X- 
ray photon indices of both XRBs and AGNs are anti-correlated 
with the Eddington ratios when the Eddington ratio is less than 
a critical value, while they become positively correl ated when 
the Eddington ratio is h igher than the critical value JWu & Gul 
120081; IGu & Caoll2008l and references therein). These results 
provide evidence for the accretion mode transition near the crit- 
ical Eddington ratio. 

It is widely believed t hat the radio emission of both the 
low/hard state XRBs (e.g. Fenderll2006l and r eferences therein 
and the low luminosity AGNs (LLAGNs, e.g.. lWu & Caol2005l 
comes from the jets. However, the detailed physical mechanism 
for the jet formation is still unclear for either AGNs or XRBs. 
The currently most favored jet for mation mechanisms include 
the Blandford-Znajek (BZ) process (Blandford & Znaiek" 19771) 
and the Blandford-Payne (BP) process (Blandford & Payne 



|1982|) . In the BZ process, energy and angular momentum are 
extracted from a rotating black hole and transferred to a remote 
astrophysical load by open magnetic field lines. In the BP pro- 
cess, the magnetic fields threading the disk extract energy from 
the rotation of the accretion disk itself to power the jet/outflow. 
Meier (199 9) showed that even if the BZ process is neglected 
entirely, the jet power contributed by the field threading the disk 
alone will still be a function of the black hole spin since the ro- 
tating metric contributes to the rotation of the magnetic field. 
Recent magnetohydrodynamic(MHD) simulations also showed 
that both the BH spin and the accretion rate in the underly- 
ing accretion disk play important roles for the j et formation 
("Koide et al.' '2000^, 'McKinnev & Gammie' '2004'; 'Hirose elalj 
2004; De VilHers et al. 2005; Hawley & KroHk 2006). The rel- 
ative importance of these two mecha nisms was e xplored by 
many authors (e.g., Gho sh & Abramowicz 11997 ; Livio et al. 



many autno rs (e.g.. luno sn & Abramowic^liyy/ lovio et al 
1999; Mei eB200UlCaol l2002; Nemmen ~l.i2007l; IWang et al 
2008). Assuming the poloidal magnetic field component at the 
disk surface to be of the same order as the toroidal field compo- 
nent, the maximal jet power extracted from the accretion disk 
(BP process) may domina te over the maximal power extracted 
by the BZ process (e.g.. lLivioetal.| [T999; Meier 2001; Caoj 
120021; lNemmenetal.1120071) . However. [Reynolds et all (120061) 
found that the dynamics of the accretion disk in plunging re- 
gion within the innermost stable circular orbit can greatly en- 
hance the trapping of large scale magnetic field on the BH, and 
therefore increase the importance of the BZ mechanism effect 
comp ared to previous estimates that ignore the plunge region 
(e.g. .Ghosh & Abramowicz. 1997.) . 

Radio galaxies are usually classified as FR I or FR II sources 
depending on their radio morphology. FR I radio galaxies (de- 
fined by edge-darkened radio structure) have lower radio power 
than FR II galaxies (defined by edge-brightened radio struc- 
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ture due to compact jet terminating hot spots; Fanaroff & Ri- 
ley 1974). What causes the morphological difference between 
FR I and FR II radio galaxies is still unclear. The theoret- 
ical models fall into two different groups: (1) the morpho- 
logical differences arise of t he different physical conditi ons in 
their ambient medium (see iGopal- Krishna & Wiital 1200ft, for 
a summary); (2) their intrinsic difference of their central en- 
gines, i.e., different accretion modes and/or jet formation pro- 
cesses (e.g., Bicknell 1995; Reynolds et al. 1996a; Meier 1999c 
Ghisellini & Celottill200U iMarchesini et al J 120041; IWold et al.) 
2007c iHardcastle et alJl2007l) ! 



Most previous works on the jet power extracted from 
ADAf s were based on self-s i milar s o lution s of ADAFs (e.g., 
Armitag e & Nat araianI 119991: iMeieiJ 120011; ICao & Rawlingsl 



20041: [Nemmen et al.ll2007h . The self-similar solution can re- 



produce the global solution quite well at large radii, while it 
devia tes significantly near the black hole (e.g., Narayan et al. 
Il997l) . where the relativistic jets are supposed to be formed. In 
this paper we calculate the jet power, incorporating some recent 
MHD simulation results, based on the global ADAF solutions 
surrounding Kerr black holes. We then compare our model cal- 
culations with the accretion power and jet power of XRBs and 
the radio galaxy FR I/FR II dichotomy. 

2. ACCRETION AND JET MODELS 

2.1. Global ADAF model 

We calculate the global structure of an accretion flow using 
an approach similar to that of Narayan et al. (1997). How- 
ever, the pseudo-Kerr potential for a rotating BH given by 
iMukhopadhvavl (l2002 l) is adopted in solving the equations of 
the accretion flow, which allows us to calculate the structure 
of an accretion flow surrounding either a spinning or a nonspin- 
ning black hole. The simple a-viscosity (t,-0 = ap, and p is total 
pressure, i.e., gas pressure plus magnetic pressure) is adopted 
and all radiative processes (synchrotron, bremsstrahlung and 
Compton scattering) are included consistently in our calcula- 
tions for ADAF structure. The advection by ions and electrons 
has been considered in the energy equation, and a more realistic 
state of accreting gas (instead of a polytropic index 7g) is em- 
ployed in the calculations, which is similar to that used by Man- 
moto (2000). We solve a set of hydrodynamical equations (i.e., 
the radial momentum, angular momentum, and energy equa- 
tions) for an ADAF, and tune the parameter k^, the specific an- 
gular momentum of the gas swallowed by the black hole, to 
let the solution passing s moothly throug h the sonic point near 
the black hole (see .Naravan et al.lll997l for details). We find 
that the derived global solutions for Kerr black holes can re- 
produce all the essential properti es of the solutions derived in 
full general relativistic frame bv iManrno to (2000) with error 
less than 10%. The global structure of an ADAF surrounding 
a BH spinning at rate j with mass Mbh can be calculated with 
proper outer boundaries (e.g., Manmoto 2000), if the parame- 
ters m, a, j3, and 5 are specified. The parameter j = J/GMgj^c~^ 
is the dimensionless angular momentum, where J is the angu- 
lar momentum of the BH, m = M/Mem is dimensionless ac- 
cretion rate, and MEdd is the Eddington accretion rate defined 
as MEdd = 1.4 X 1O'^Mbh/M0 g s"'. The value of a adopted 
in ADAF modeling is suppose d to be within a very narrow 
range, i.e., a = 0.1-0.3 (e.g.. iNaravan & McClintockl l2008l 
and references therein), which is supported b y MHD numerical 
simul ations of accretion flows ~ 0.05-0.2 dHawlev & BalbusI 
l2002h and the observationally-determined values ~ 0.1—0.4 



based primarily on studies of outbursts in dwarf novae and 
X-ray transients) King 2007). The magnetic parameter (3 (de- 
fined as ratio of gas to magnetic pressure in the accretion flow, 
/3 = Pg/Pm) is not an independent parameter and can be related 
to a as /3 ~ (0.55 -a)/o ;, as suggested by MHD simulations 
(e.g.. lHawlev et al.lll995"l) . where P^ = Bdynamo/^T^ and Bdynamo 
is the magnetic field strength of the ADAF in the local refer- 
ence frame. The parameter /? ~ 1 -5 for the typical value of 
a ~ 0. 1 -0.3. Another poorly constrained parameter is S, which 
describes the fraction of the turbulent dissipation that directly 
heats the electrons in the flow. Recent ADAF models typically 
assume S ~ 0.3-0.5 (e.g.. Yuan et al. 2003; Wu et al. 2007, see 
also Sharma et al. 2007 for slightly lower 6 value). 

Two important modificat ions of the above global ADAF 
model are also included (see lMeienl200ll : lNemmen et al.ll2007l 
for more details). First, as viewed from an outside observer at 
infinity in the Boyer-Lindquist reference frame, the disk angu- 
lar velocity fi is a sum of its angular velocity relative to the lo- 
cal metric fl plus the angular velocity of the metric itself in the 
Boyer-Lindquist frame tu = -g(t>t/g4>4>, i.e., fi = fJ +ijj. Second, 
we also take into account the field-enhancing shear caused by 
frame dragging in the Kerr metric, as first suggested bv iMeien 
^999), w hich see ms to be supported by MHD sim ulations (e.g ., 
iHawlev & Kroli k'2b06). Foflowing the work of Meier (2001]), 
the amplified magnetic field related to the magnetic field pro- 
duced by the dynamo process in the ADAF can be expressed as 
B = gBdynamo, whcrc g = Vt/Vt is the field-enhancing factor. 

2.2. Evaluating the jet power of Blandford-Znajek model 

For a black hole of mass Mbh and dimensionless angular 
momentum j, with magnetic fields B±^ normal to the horizon 
at Rii = [l+{l-fy^^]Rg (Rg = GMbb/c^ is the gravitational 
radius), t he power extracted with the BZ me chanism is given 
by (e .g., iGhosh & Abramowicd 119971; iMacDonald & Thornd 
119821) 

Qlf=^colBlRifc, (1) 
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where Wp = riF(^H- ^f)/^h i^ determined by the angular ve- 
locity of field lines Up relative to that of the hole Hh- In or- 
der to estimate the maximal power extr acted from a spinning 
BH, u jp is always required to be 1/2 (e.g.. lLivio et al.lll999l; ICaol 
I2OO2I) . The magnetic field B± is assumed to app roximate to the 
poloidal component Bp, and iLivio et al] (1 19991) proposed that 
Bp ~ Bdynamo duc to the hot thick disk of ADAF (H r^R). There- 
fore, we use B± ~ Bp ~ ^Bdynamo in our calculations consider- 
ing the field enhancing effect. This is consistent with recent 
MHD simulations in which t he poloidal fields are dominant in 
the z-dir ection near the BH (|Kato et al.ll2004l) . Following the 
work of 'Nem men et al.l (l2007l) . all the physical quantities are 
evaluated atR = /?ms- 

2.3. Evaluating the jet power of hybrid model 

As pointed out bv lMeierl (Il999l) . the differential dragging of 
the frames will also act as a dynamo to amplify the magnetic 
field, and therefore even if the BZ process is neglected entirely, 
the jet power contributed by the field threading the disk alone 
will still be a function of the BH spin, since the metric of the 
rotating black hole contributes to the rotation of the magnetic 
field. Both the BP and BZ mechanisms were incorporated in 
this hybrid jet formation model, in which the magnetic fields 
extract energy both from the accretion flow and the spinning 



hole dMeieJl 99911200 iHPunslv & Coronitil 19901) . The total jet 
power for the hybrid model is given by 
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e]if = bX"V32c, (2) 



where R is the characteristic size of jet formation region, Bp : 

^Bd mamo, O = f} +a;, and all quantities are evaluated at R = Ra 
(see lMeieiJl200lL for details). 



3. RESULTS AND DISCUSSION 

To estimate the jet power extracted from the inner region of 
the disk, we employ the global ADAF solution for a spinning 
black hole, and the field-enha ncing shear in the Kerr metric has 
been taken into account (e.g..' Meienl2001h . Our global calcu- 
lations show that the magnetic field can be amplified roughly 2 
times when j ^ 0.9 in the plunging region (the zone between 
Rms and Ru)- We find that the jet power for either the BZ 
model or the hybrid model is roughly proportional to the ac- 
cretion rate/BH mass, but its dependence on j is rather com- 
plicated. Figure 1 shows the spin dependence of jet power for 
the BZ model and the hybrid model for given Mbh = 10^ Mq 
and m = 0.01 (also see similar results in Fig. 1 of Nemmen 
et al. 2007). The black-solid lines and blue-dashed lines in 
Fig. 1 denote the jet power for two different values of the vis- 
cosity parameter a = 0.3 (/? ~ 1) and a = 0.1 (/3 ~ 5) for the 
case of (5 = 0.5, respectively. We find that the jet power varies 
little for different viscosity parameters, provided all other pa- 
rameters are fixed. The physical reason is that the jet power 
(BZ/hybrid models) Qjet oc B^ oc 1/q!(1 +/3) ~ constant in the 
ADAF (see Narayan & Yi 1995, and relation of a and /3 in Sect. 
2.1). The red dotted lines denote the jet power for the cases of 
a = 0.3 (/3 ~ 1, Fig. 1) and (5 = 0.1, which indicates that the 
jet power is also not sensitive to the value of S. We find that 
the ratio of the gas pressure to the magnetic pressure ~ 0.4- 1 
in the plunging region varying with different BH spin param- 
eters j is mainly due to the field-enhancing shear caused by 
frame-dragging. Our results for the case of a = 0.3 are roughly 
consistent with the MHD simulations in Kerr metric indicating 
the ratio ^ 0;3 - 1 at the inner boundary (iDe Villiers et al.l2003l: 
iHirose et al.ll2004| ). Therefore, a = 0.3 and S = 0.5 are adopted 
in the following calculations. We find that the values of these 
parameters will not affect our main conclusions. 

We find that the jet power of the hybrid model, Qf^f^, is nearly 

25 times larger than that of the BZ model, Qf^, for a rapidly 
spinning black hole with j ^ 0.9, or even higher for a smaller 
j (Fig. 1). The jet efficiency [defined as TTjet = 2jet/A/(/?ms)c^] 
for a spinning BH with j = 0.99 is 0.8% and 20% for the BZ 
model and the hybrid model, respectively. This suggests that 
the hybrid model plays a more important role than the pure 
BZ model in the jet formation for a black hole surrounded 
by an ADAF. Our results of the hybrid model are consistent 
with 77ipf ~ 21% for / = .99 in the numerical MHD simulations 
dHawlev & Krolikll2006l) . Hereafter, we implicitly consider the 
jet power Qjet — Qtf^ in our calculations, since the jet power 
of the BZ model is always much lower than that of the hybrid 
model for ADAFs. 



3.1. Accretion Power/Jet Power of ADAF and Comparison 
with XRBs 
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Fig. 1 . — The relation between the jet power Qjet extracted from the under- 
lying ADAF and BH spin j for the BZ and hybrid models, respectively. The 
black solid lines and blue dashed lines denote a = 0.3 (/3 = 1) and a = 0. 1 
(/3 = 5 and 5 = 0.5), respectively. The red dotted lines denote a = 0.3 for the 
case of (5 = 0.1. The BH mass Mbh = 10*Mq and accretion rate m = 0.01 are 
adopted in the calculations. 

Once accretion rate m falls below a critical value ma-, the 
standard disk will transit to an ADAF at radii less than the tran- 
sition radius. In the ADAF scenario, the accretion power (or 
radiated luminosity) is much lower than that of SSD as a re- 
sult of reduced radiative efficiency. The power created by ac- 
creting matter can either be radiated away, advected into the 
BH, or taken away by jet/wind. Normally, the radiative effi- 
ciency of the pure ADAF is a function of accretion rate, i.e., 
77 = rjim/rhti)^. The radiative efficiency rj of the ADAF is 
roughly proportion to the accretion rate (e.g., ^ ^ 1) at lower ac- 
cretion rate and sligh tly flattens (^ < 1 ) at ne a r the critical value 
rhtr JM erloni et al. 2001 lYuan & Cuil I200I IWu & Cad I2OO6I: 
ISharma et alji2007,) . However, the radiative efficiency should 
be slightly higher when the accretion rate is near m^, especially 
considering the contribution of the outer SSD which should be 
important when the transition radius is s mall, and/or the possi- 
bility of a condensation-feed inner disk dLiu et alj|2007h . Here, 
we do not consider these two mechanisms in details due to both 
of them are still quite unclear We assume that ^ = 1 when 
m ^ nia (see also Narayan & McClintock 2008), whi ch is con- 
sisten t with that constrained from the XRBs (e.g., Kording et al.l 
I2OO6I) . and will not affect our conclusion. Therefore, the bolo- 
metric luminosity (or accretion power) of the ADAF can be de- 
scribed by 



OTtr 



(3) 



where r](j) is radiative efficiency depending on the BH spin 
parameter j. The exact value of th e critical accretion rate 
witi is still unclear. iMaccarong (l2003l) found mtr '^ 0.02 from 
a detailed investigation of state transitions in X-ray binaries. 
IWu & Gul (|2008.) suggested that the accretion rate for the disk 
transition may be 2-3 times lower than that derived from the 
state transition. We adopt mtr = 0.01 in our calculations. Fig- 
ure 2 shows the relation of 2fet''/^Edd - ni (solid line) and 
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Lboi/LEdd-rn (dashed line) for j = 0.99,0.9,0.7, and 0, respec- 
tively. We find that the jet power dominates over the accre- 
tion power when the accretion rate is less than a critical value 
m < nic and the objects enter into the "jet-powe r-dominated 
advective systems," where 3 x 10"^ ^ "Jc ^ 5 x lO""* depends 
on the BH spin parameter j. The accretion power is dominant 
when m > nic, and the objects are therefore ''accretion-power- 
dominated advective systems." We note that the uncertainties of 
5 parameter will lead to slightly different critical accretion rate 
rhc due to the radiative efficiency of ADAF, even the jet power 
is not sensitive to the parameter 5. We find the uncertainties 
on the critical accretion rate mc will be less than 5 times when 
considering ^ ~ 1 ±0.2 for 0.1 < 5 < 0.5 and BH spin ; = 0- 1, 
and will less th an 2 times for th e typical value j = 0.7- 1 found 
in XRBs (e.g.. iLiu et alj|2008l and references therein), which 
will not affect our main conclusion. 
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Fig. 2. — The relation between Qjct/iEdd" 
(solid line) for BH spin parameters: j = 0.9 
respectively. 

Our calculation show that the jet power 2jet oc «' oc m 
and accretion power Lboi oc m^. The relation gjet oc L!^-^^ is 
consistent with the nonlinear correlation between radio lumi- 
nosity and X-ray luminosity, Lradio oc L^^, in the low/hard 
state (e.g.. ICorbel et al.l 120031) . when considering the simple 
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optically thick conical jet model, where Qjm oc L i„ (e-g-, 
iBlandford & Konigllll979HFalcke & Biermaiinlll995 b and Lx is 



an indicator of Lboi- iFender et al. 



(|2003|) gave a conservative 
lower limit of the critical accretion rate m^ ^ 10""* for the tran- 
sition of "accretion-power-dominated regime" and "jet-power- 
dominated regime" assuming the power created by the accreting 
matter is radiated (accretion power) and tak en away by jet (jet 
power ) based on the XRB XTE Jill 8-H480. rMighari & FendeJ 
(|2006|) further suggested that these two regime transitions may 
occur at the slightly higher critical value rric ^ 10"^, using a 
slightly higher jet power normalization. These estimates based 
on low/hard state XRBs are roughly consistent with our re- 
sult that the critical accretion rate is 3 x 10""* < m^ ^ 5 x lO"-' 
when considering the uncertainties in the estimates of the jet 
power normalization from the observation. We should note that 
the energy advection play an important role for jet formation. 



since only a small fraction of the gravitational energy of the 
flow is dissipated locally, and most of the energy heats the pro- 
tons/electrons and is advected inward. This leads to a hot thick 
disk, which allows high poloidal magnetic field strength in the 
inner region of the flow. We find that, as an example, the frac- 
tion of output energy ^ (gjet+LfioO/O.lMc^ '^0.2 for m= lO"-' 
and j = 0.7, and the other 80% of the energy is advected into 
the central BH. Therefore, the BH central engine will be a 
"jef-/?ower-dominated advective system" when the luminosity 
(or accretion rate) is less than a critical value, which is con- 
sistent with that constrained with the observations on the BH 
XRBs and neutron stars (e.g., Kdrding et al. 2006). It should be 
noted that we use "jet-/?ower-dominated" not "jet-dominated" 
to discriminate the possible confusion about the origin of multi- 
wavelength emission. For example, the X-ray emission in the 
low/hard state XRBs and LLAGNs is still controversial, and 
it may be dominat ed by the jet emission (e.g., Markoff e t al.l 
120031; iFalcke et all |2 004). by flie underlying ADAF, or both 
(e.g.. lYuan & Cuill2005.:.Wu et^l2007h . 

3.2. FR I/II Dichotomy 

The division between FR I and FR II radio galaxies is clearly 
shown by a line in the plane of the total radio luminosity and th e 
optical luminosity of the host galaxy ( Ledlow & Owenlll996l) . 
iGhiselUni & Celottil (1200 1.) used the optical luminosity of the 
host galaxy and the radio luminosity to estimate the mass of its 
central BH and jet power, respectively. They proposed that the 
FR I/FR II separation can be interpreted by the systematically 
different ratios of the jet power to BH mass for FR I and FR II 
sources. This implies that the FR I/FR II division is linked to the 
physics (accretion and/or jet processes) on very small scales. 
They argued that if the jet power is related with accretion power, 
the accretion mode for low-power FR I sources may be dif- 
ferent from that for powerful FR II sources. There is grow- 
ing evidence to suggest that most FR I type radio galaxy nu- 
clei may possess AD AFs (Reynolds et al. 1996b; Gliozzi e tal.l 



2003; Merloni et al. 2003; Do nato et al. 2004; Wu et al. 2007). 
Cao & Rawlings (2004) found that the jet power extracted from 



the underlying ADAFs with pure BZ mechanism is insufficient 
for some high power-FR I radio galaxies. In this work, we com- 
pare the jet power of the hybrid model with the observed FR 
I/FR II dichotomy in the Mbh - fijet plane. 

The dividing li ne be tween FR I and FR II sour ces of 
Ledl ow & OwenI d 19961) is given approximately by iMeierl 
(il99C as 

logP,ad = -0.66MR+ 10.35, (4) 

where f,ad is the observed radio power at 1.4GHz (in W Hz"') 
and Mr is the absolute optical /?-band magnitude of the host 
galaxy. An empirical relation between Mr of the host galaxy 
and central BH mass, 

log(MBH/M0) = -0.50(±0.02)M^-2.96(±0.48), (5) 

was derived by McLure & Dunlop (2002). The jet power is usu- 
ally estimated from the radio luminosity by using the relation 

gjet - 3 X \Q'^fl^L%\ W, (6) 

where L151 is the total radio luminosity at 151 MHz in 
units of lO^'^W Hz"' sr"' and the factor / parameterizes 
the uncertainties of the nor malization which is constrained 
to be between 1 and 20 (se e IWiilott et al.|[T999l for details). 
iBlundell & Rawlingsl (1200ff) argued / to be most likely in the 
range of 10-20. The radio luminosity at 1 .4 GHz is converted to 
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the luminosity at 151 MHz by assuming a radio spectra index 
a = 0.8, where L{i^) oc i^~". Thus, we can obtain the dividing 
line between jet power and BH mass by using the Eqs. (4)-(6), 

logGje,(ergs-i)= 1.131ogMBH(M0) + 33.42+1.5Olog/. (7) 
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Fig. 3. — Jet power and BH mass divid ing line between FR I and FR II radio 
galaxies for the ILedlow & Owen j 19961) sample (solid lines; the jet power is 
calculated from the 151 MHz radio luminosity with two different correction 
factors / = 20 (top) and / = 10 (bottom), respectively). The dashed and dotted 
lines are the jet power of the hybrid model and BZ model extracted from the 
underlying ADAFs with j = 0.99,0.9,0.7 {from top to bottom), respectively. 
The accretion rate m = 0.01 is adopted in all calculations, which is roughly 
consistent with that constrained from the ionization luminosity of the dividing 
line between FR I and FR II sources (see the text for more details). 

The maximal jet power of the hybrid model {solid line) and 
the BZ model {dotted line) for ADAFs are plotted in Fig. 3 
with different black hole spin parameters: j = 0.99,0.9,0.7, re- 
spectively. The accretion rate m = »% =0.01 is adopted in our 
calculations. We find that the maximal power of the pure BZ 
model is nearly 25 times less than the dividing line between FR 
I and FR II radio galaxies. However, we find that the maximal 
jet power of the hybrid model with j = 0.9/0.99 can roughly re- 
produce the jet power dividing line with / = 10/20 (see Fig. 3). 
Thus, the dividing line between FR I and FR II sources corre- 
sponds to the maximal jet power of the hybrid model extracted 
from the ADAFs surrounding rapidly spinning black holes ac- 
creting at the critical rate riia = 0.01, which is consistent with the 
analogue between FR I sources and the low/hard state of XRBs 
(e.g., Merloni et al. 2003; Falcke et al. 2004). We find that the 
jet power of the hybrid model is very sensitive to the BH spin, 
especially when j > 0.9 (Fig. 1). Therefore, it seems that this 
sensitivity would likely blur the FRJ/FR II dividing line if the 
spins of black holes spread around 0.9. The FR I/FR II dividing 
line is indeed not very clear (one can find a few FR I sources 
above the dividing line , and vice versa for FRII sources, e.g., 
iLedlow & OwenI (11996 !)). The critical accretion rate mtr = 0.01 
for the accretion mode transition is also supported by the ion- 
ization luminosity for t he sep aration between FR I and FR II 
sources. IWillott et al] (Il999h found that the ionization lumi- 
nosity of radio galaxies is roughly equal to the jet power for 



/ = 20, which corresponds to Lion/i'Edd ^ 2. 5 x 10"^ for typi- 
cal BH mass Mbh ~ lO^^'^^M© in sample of lLedlow & OwenI 
(Il99d) . This ionization luminosity Eddington ratio roughly cor- 
responds to m ~ O.OlQ if we assume that the bolometric lumi- 
nosity is equal to the ionization luminosity and the BH is spin- 
ning rapidly [e.g., ri{j) > 0.2 for j > 0.95]. It is still unclear 
why the jet power of FR II radio sources is always above this 
dividing line, which is beyond the scope of this work. 

Our calculations on the jet power for either the BZ or BP 
mechanisms are based on the pure ADAF model. iNaravan & Yil 
d 1 9941) found that the ADAF has a positive Bernoulli param- 
eter, and the accretion flow is therefore gravitationally un- 
bound, which implies that the gas may escape as outflows. 
In this case, the accretion rate of the advection-dominated 
inflow-outflow solution (ADIOS) is a function of the radius 
instead of a constant accretion rate for the pure ADAF (e.g., 
Iglandford & Begelman 1999). For ADIOS, the gas swaflowed 
by the BH is only a fraction of the rate at which it is sup- 
plied, as part of the gas is carried away in the winds before 
it reaches the BH. Therefore, both the accretion power and jet 
power will decrease in the ADIOS compared the pure ADAF 
case for given accretion rate at the outer boundary. The accre- 
tion power Lboi is roughly oc ;«?,, while the jet power Qjet is 
roughly oc niin at the inner edge of the flow. It means that, in 
Fig. 2, the lines for Lboi shift down more than those for Qjet 
in the presence of winds. Thus, the critical accretion rate nic 
for the transition between ''jet-powe r-dominated advective sys- 
tems" and "accretion-power-dominated advective systems" in 
the ADIOS should be slightly higher than that derived from the 
pure ADAF case. Therefore, the critical accretion rate in the 
ADIOS maybe close to m^ ~ 10"^ for the case of j > 0.9, which 
is roughly consistent with that constrained from the XRBs (e.g., 
iMigliari & Fenderfl2006l) . The ADIOS has a structure and a 
upper limit on the accretion rate at the inner edge of the disk 
simila r to those of a pure ADAF if the wind is not very strong 
(e.g., Chang et al.l2002l:lNemmen et al.l2007h . The jet power of 
BZ/BP mechanisms is dominantly extracted from the inner re- 
gion of the flow, so the maximal jet power for an ADIOS should 
be similar to that for an ADAF without winds, provided that 
their accretion rate at the inner edge of the flow are the same. 
Therefore, the maximal jet power from ADIOS with mi,, ~ 0.01 
and /■ ^ 0.9-0.99 can still reproduce the dividing line of the FR 
I/FR II dichotomy (Fig. 3). 

4. SUMMARY AND CONCLUSION 

The main conclusions of this work can be summarized as 
follows: 

(1) We investigate the jet power of the BZ mechanism and 
hybrid mechanism for the ADAFs surrounding rotating black 
holes based on our global ADAF solutions in Kerr metric. 
We find that the jet power of the hybrid model is about I 
order of magnitude higher than that of the BZ model, and 
the jet efficiency of the hybrid model is roughly consistent 
with that of the numerical MHD simulations for Kerr black 
holes(Hawley & Krolik 2006). 

(2) The jet power dominates over the accretion power when 
the accretion rate is less than a critical value m^, while the 
accretion power will be dominant when the accretion rate is 
larger than this critical value (Fig. 2), which is roughly consis- 
tent with that constrained from the observations of XRBs (e.g.. 



* We note that the differences of the Eddington ratio in Ghisellini & Celotti (2001') and lWold et all (2007) from ours are caused by using different factor /, and / = 1 
is adopted in their work, and some coefficients of Eqs. iniGhiselhni & Celotti, t200U) may be erroneous as pointed out bv iWold et alj 1 ,2001,) . 
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iFender et al.ll2003MMigliari & Fendeij|2006h . 

(3) The dividing line of the Ledlow-Owen relation in Qjet - 
Mbh plane for the FR I/FR II dichotomy can be well reproduced 
by the hybrid jet model, provided that the ADAFs surrounding 
Kerr black holes are accreting at m < 0.01 in FR I sources(Fig. 
3). The ionization luminosity corresponding to the dividing line 
plane also suggests the critical accretion rate to be ~ 0.01 if the 
empirical relation between the photo-ionizing luminosity and 
the radio luminosity is adopted. These strongly imply that most 



FR I radio galaxies may have a different accretion mode from 
FR II sources. 
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